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Dynamical mean field methods are used to calculate the phase diagram, many-body density of 
states, relative orbital occupancy and Fermi surface shape for a realistic model of LaiViOa-based 
superlattices. The model is derived from density functional band calculations and includes oxygen 
orbitals. The combination of the on-site Hunds interaction and charge-transfer between the transi- 
tion metal and the oxygen orbitals is found to reduce the orbital polarization far below the levels 
predicted either by band structure calculations or by many-body analyses of Hubbard-type models 
which do not explicitly include the oxygen orbitals. The findings indicate that heterostructuring 
is unlikely to produce one band model physics and demonstrate the fundamental inadequacy of 
modeling the physics of late transition metal oxides with Hubbard-like models. 

PACS numbers: 73.21. Cd, 73.20.-r, 71.10.-w 



The electronic properties of transition metal oxides are 
of central importance to condensed matter physics and 
materials science, both for their fundamental scientific 
intereslP^] and for their potential for novel applicationsP 
Of particular current interest are the new possibili- 
ties enabled by advances in atomic-scale layer-by-l ayer 
growth of combinations of complex oxide materials^. 
Experiments report remarkable and unexpected prop- 
erties including interface superconductivity^ orbital 
reconstruction^ high Curie-t emp erature magnetisuP 
and metal-insulator transitions 

These developments suggest that it may become pos- 
sible to design materials with desired 'correlated elec- 
tron' properties such as high-temperature superconduc- 
tivity. Rational materials design requires knowledge of 
the structure-property relation between atomic arrange- 
ment and electronic density of states. In systems such as 
conventional semiconductors where the electronic struc- 
ture is well described by density functional band theory 
this issue is well understood, so that e.g. band-gap engi- 
neering by choice of superlattice or size of quantum dot 
is now routine.^ However, band theory provides an in- 
complete description of the relevant electronic states in 
materials with strong electronic correlations, and as a re- 
sult the understanding of structure-property relations in 
transition metal oxides is much less well developed. 

The important states transition metal oxides are the 
transition metal d-orbitals. In the CtiC^-based high- 
ly cuprates the nominal electronic configuration of the 
Cu is d 9 (one hole in the d-shell) and the crystal struc- 
ture breaks the rotational symmetry to the point that 
the only relevant Cu orbital is the d x 2_ y 2. The result- 
ing quasi-two-dimensional x 2 — ?/ 2 -dominated electronic 
structure is belie ved t o be crucial to the high transi- 
tion temperature! 12 * 13 -^ In the related material LaNiOs 
the nominal configuration is Ni d and the relevant or- 
bitals are d x 2_ y 2 and d^ z i_ r i^ which transform as a dou- 
blet under crystal symmetry operations. A recent pa- 



per suggested that in a superlattice composed of alter- 
nating layers of LaNiO^ and a related but insulating 
material such as LaAlOj, the d 3z 2_ r 2 orbital would be 
pushed so far away in energy that the Fermi surface would 
only have one sheet possibly leading to high tempera- 
ture superconductivity.^ While there is to our knowl- 
edge no experimental evidence for superconductivity in 
nickelate superlattices, essential aspects of correlated- 
electron behavior including metal-insulator transitions 
and magneti c be havior are intimately connected to or- 
bital physicspES so th e i ssue f orbital polarization in 
superlattices and its relation to electronic behavior is an 
essential question for the subject of 'correlated electron 
materials by design'. _™ 

Density functional band theory calculations^^ indi- 
cate that an LaNiO^/ LaXOy, superlattice has partial 
(~ 20 — 40%) polarization of the Ni d orbitals, with 
the precise amount (and sign) of the polarization de- 
pending on the specific choice of counter-ion X in the 
heterostructurep^ but the known inadequacy of band 
theory for charge transfer materials such as the rare earth 
nickel oxides motivates an examination of the effects of 
correlations. In this paper we show that strong correla- 
tion effects actually decrease the polarization. The key 
feature aspect of our calculation is the use of a realis- 
tic Hamiltonian which is derived from the density func- 
tional band calculations of Ref.[T7|and in particular takes 
oxygen states into account explicitly. The possibility of 
charge transfer to the oxygen affects the results in an 
important way. The model is of the general form 

H = Hd + Hhyb + Hu gan d (1) 

The Hamiltonian for the correlated ("d") subspace is 
H *=Y, ^d) aa d 3aa + |iVi (Ni l) + Hj (2) 

with a = 1, 2 labelling the x 2 — y 2 and 3z 2 — r 2 d states, 



2 




A(eV) 



FIG. 1: Metal-insulator phase diagram of two dimensional 
Ni — O superlattice calculated using single-site dynamical 
mean field theory in plane of on-site interaction strength U 
and charge transfer energy A = e p — Ed- Plotted is U c i, the 
lower limit of stability of the insulating phase 



Ni the total number of d electrons on site j, U the on- 
site repulsion and Hj the additional 'Hunds' interactions 
which give the multiplet structure at given d occupa- 
tion. Consistent with recent downfolding studies of the 
screened interacti on m atrix elements in transition metal 
oxide compounds! 18 * 19 ^ we take the interaction terms to 
have the standard Slater-Kanamori formPset J = 0.5eV 
and consider a range of U. 

Hhyb and Hu gan d describe the d — p hybridization and 
the ligand (oxygen) part of the band structure. They 
may be obtained by fitting a tight binding model to a 
calculated band structure or, equivalently, by construct- 
ing Wannier functions from st ates within some appro- 
priately chosen energy window? 20 * 21 ! The resulting tight 
binding model (obtained from non-spin-polarized GGA 
band calculations^) involves e g orbitals as well as oxy- 
gen 2p,j orbitals with a mean oxygen energy e p , oxygen- 
oxygen hoppings t pp and Ni-0 hoppings t p d and is de- 
scribed in detail in Ref. |T7] The choice X = In was found 
to produce the largest d x 2_ y 2 occupancy and a very weak 
hybridization in the direction transverse to the superlat- 
tice plane. In this paper we study the extreme case of 
a two dimensional system defined using the Ni and O 
parameters obtained from the LaNiO^/ LalnO^ super- 
lattice studied in Ref. [17] but with the hopping from the 
apical O ion to the X plane set to zero. 

The charge transfer energy A = e p — Ed plays a cru- 
cial role in the theory of the metal-insulator transition in 
oxides! 2 ^ A is renormalized from the band theory value 
by the double-counting correction^ whose value is an 
important unsolved problem in mate rials th eory. Differ- 
ent prescriptions have been proposed) 20 * 21 * 23 ! but no clear 
consensus has emerged. We study a range of A here, 
corresponding to a range of double counting corrections. 

To treat the many-body physics we use the single- 
site dynamical mean field approxima tion^ with the 
hybridization expansion impurity solveJ^ISll which can 



treat the full rotationally invariant Hunds coupling. We 
focus mainly on the metallic regions of the phase dia- 
gram. Care must be taken to converge the solution (up 
to 30 iterations of the dynamical mean field procedure are 
required for parameters near the metal-insulator transi- 
tion line) and temperatures must be chosen low enough 
to reveal the quasiparticle behavior in all orbital sectors. 

In the single-site dynamical mean field approximation 
the metal-insulator transition is first order, character- 
ized by an upper critical interaction strength U C 2 which 
marks the limit of stability of the metallic phase and a 
lower critical interaction strength U c i which marks the 
limit of stability of the insulating phase. Fig. [I] presents 
U c i of A. The location of the rare earth nickelates on the 
phase diagram is not known. LaNiO^ is metallic in bulk 
and (except for one and perhaps 2 monolayer samples) 
in thin film form. Other members of the family such as 
NdNiOs have insulating ground states, suggesting that 
the materials are close to a metal/charge-transfer insu- 
lator phase boundary, but the origin of the insulating 
phase is controversialF^nSSl We therefore study a range 
of parameters in the metallic state. In DMFT the insu- 
lating state typically exhibits some form of orbital order, 
making the interpretation less clear. 

The many-body electronic structure is represented by 
the local spectral function (many-body density of states) 
A a {uj) = Im(-iJ dte-^ ([Mj,t)^lU,t')])) (here a 
labels an orbital and j a unit cell in the lattice). Panel 
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FIG. 2: (Color online) Analytically continued spectral func- 
tion A(uj), for Ni d x 2_ y 2 (dotted, blue on line), d 3z 2_ r 2 (solid, 
red on-line) orbitals and sum of all O orbitals (dashed, ma- 
genta on line). Parameters (Nd computed from integral of the 
many-body DOS over the entire manifold of occupied states) 
(a) U = 0, N d = 1.99. e d = -1.22 eV, e p = -5.2 eV. (b) 
U = 4 eV, N d = 1.60. e d = -3.91 eV, e p = -7.89 eV. (c) 
U = 6 eV, N d = 2.12. s d = -8.95 eV, e p = -4.93 eV. (d) 
U = 6 eV, N d = 1.52. e d = -5.75 eV, e p = -7.73 eV. 
Fermi level is zero. J = 0.5 eV and T = 0.1 eV. Com- 
puted polarization P from Eq. [4] with Ei ow — — 3eV 
shown on each panel. The P values corresponding to inte- 
gration over the energy range of the entire p — d manifold 
areP = (0.16,0.14,0.11,0.11) for panels a-d respectively. 
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FIG. 3: Calculated Fermi surfaces at J — 0.5 eV and A = 
-1.98 eV. (a) U = 0, N d = 2.44, e d = -1.81 eV. (b) U = 4 
eV, N d = 1.65, e d = -4.65 eV. (c) U = 6 eV, 7V d = 1.45, 
e d = -5.78 eV. (d) U = 7 eV, iV d = 1.34, e d = -6.02 eV. For 
panels (a)-(c) T = 0.05 eV; for panel (d) T = 0.025 eV, 



(a) of Fig.[2]presents the noninteracting (U = 0) case, for 
which A w —4eV. These spectra are consistent with pre- 
viously published GGA results^ (the small differences 
arise from the difference between the two dimensional 
model used here and the fully three dimensional model 
of Ref. I17j) . Two energy regions are evident: a near 
Fermi surface region representing the d — p antibonding 
band and a lower energy region representing the bonding 
combination of d and p. The p level energy is visible as 
a sharp peak in this energy region. Panel (b) shows the 
result of increasing the interaction strength to U = AeV 
while keeping both the total electron count and A fixed. 
Comparison to Fig. [I] shows that this change moves the 
system close to the metal-insulator phase boundary. We 
see that the splitting between the bonding and antibond- 
ing regions of the spectrum increases, essentially because 
the interaction increases the effective d-level energy. We 
also see that the differences between the spectra of the 
two d orbitals are smaller than in panel (a). Panels (c) 
and (d) of Fig. [2] shows spectra obtained for a stronger 
interaction U — 6eV. Panel (c) shows a double-counting 
correction corresponding to A = 4eV and chosen to undo 
the shift in the d-p splitting. The parameters are far 
from the metal-insulator phase boundary and the spec- 
tra are seen to be very similar to those computed for the 
noninteracting model. In panel (d) we have chosen the 
double counting correction to keep the energy separation 
between the O states and antibonding p—d states approx- 
imately the same as in panel (b). The features near the 
Fermi level are narrower than in panel (b) because the 
system is closer to the metal-insulator phase boundary, 
but the spectra are otherwise similar. 



The occupancy n a for a given orbital a is defined as 

dcu 

n a = / 

JE, ow 7T 

and the orbital polarization is 



-A a {oj) 



(3) 



P = 



n x 2_ y 2 — n 3z 2_ 

11x 2 -y 2 + n 3z 2 ~ 



(4) 



We have chosen the zero of energy such that \x = 0. We 
believe it is most physically reasonable to focus on the 
difference in occupancy of near Fermi-surface states, cor- 
responding to taking Ei ow = —3eV to capture the anti- 
bonding but not the bonding bands. Alternatively, one 
may integrate over the whole (many-body) bandwidth. 
We have provided the P corresponding to both defini- 
tions; the results are very similar and the conclusions are 
not changed. The computed d occupancies and polar- 
izations are given in the caption and panels of Fig. [2j 
interactions decrease the polarizations. 

An alternative, low-energy definition of orbital polar- 
ization may be obtained from the Fermi surface. In pseu- 
docubic LaNiO^ the calculated Fermi surface has two 
sheets, corresponding to the two relevant d orbitalsP^ 
In single-layer high-T c cuprates the Fermi surface has 
only one sheet, corresponding to a single relevant d band, 
so one may identify a single-sheeted Fermi surface with 
an orbitally polarized low energy theory. Fig. [3] shows 
the evolution of the Fermi surface of the nickelate super- 
lattice as the interaction strength is increased at fixed 
A « -2. The non-interacting model has a substantial 
degree of orbital polarization, as seen from the very small 
size of the central Fermi surface region, but as soon as the 
interaction is turned on the size of the central patch in- 
creases and then does not change over the entire metallic 
region, consistent with the values of P given in the cap- 
tion of Fig. [5] (The Fermi surface in panel (d) is slightly 
smaller because at the lowest accessible temperature the 
fully coherent Fermi liquid state was not achieved.) 

The small value of P we find is in agreement with re- 
cent resonant X-ray absorption experiments^ but does 
not agree with results of previous dyn amica l mean field 
studies of Hansmann and collaborators While there 
are minor technical differences (including the use, by 
Hansmann ct. al. of an Ising approximation to the 
Hunds interaction) we believe that the most important 
issue is the model. Refs. H6|32l downfolded the band the- 
ory results to a two-band model representing only the 
antibonding band, whereas in our work the Ni-O charge 
transfer plays an important dynamical role, enabling the 
high spin d s L configuration which is not susceptible to 
orbital polarization. 

Our results suggest that in realistic models of nicke- 
late superlattices, a significant orbital polarization will 
be very difficult to achieve. However, the two orbitals 
will not have identical properties. We present in Fig. [4] 
the difference AA(w) = A x 2_ y 2(uj) — A 3z 2_ r 2(uj) calcu- 
lated for parameters corresponding to Fig. [2]ja)-(d). The 
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FIG. 4: (Color online) d-spectrum difference plot AA(oj) = 
A x 2_ y i(u)) — A 3z 2_ r 2(uj) (solid line, red online). Parameters 
for (a)-(d) correspond to those in Fig. [2] 

d-spectral function may be measured in resonant x-ray 
scattering experiments, and difference spectra are rela- 
tively insensitive to experimental complications such as 
core-exciton and final-state corrections. We see that the 
different electronic structures lead to observable effects 
on the spectra. In the ~ 2.0 case (corresponding to 
panels (a) and (c) in Fig. [2]) the difference spectra reveal 
prominent peaks just below the Fermi level (~ —0.5 eV) 
and around ~ +2.0 eV. These features do not appear in 



the Nd ~ 1-45 cases, (panels (b) and (d) in Fig]2b. Also 
the two feature provide a measure of the physical e p — e ( i- 
Comparison of these calculat ions to new generations of 
X-ray absorption experiments^^ may help evaluate the 
orbital polarization and pin other material parameters. 

In summary, we have shown that in a realistic many- 
body model of nickelate heterostructures, it is essentially 
not possible to achieve a significant degree of orbital 
polarization, so that the idea^ of obtaining a single- 
band electronic structure must be discarded. Further, 
we showed that a r educt ion of the full Hamiltonian to a 
Hubbard-like modePSES which includes only the corre- 
lated orbitals yields a fundamentally misleading picture 
of the electronic structure. We presented spectra which 
should help in establishing the actual value of the double- 
counting correction for these materials, which is crucial 
to the metal-insulator transition behavior. 
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